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1 Introduction 



In the electroweak sector of the Standard Model (SM) every particle acquires its mass 
through an interaction with a scalar potential in a non-trivial vacuum. As a consequence, 

— 1/2 

all the masses are proportional to a common scale, namely G F , which is fixed by low- 
energy measurements such as the p, decay rate. In this situation, the decoupling theorem 
[T] does not hold and thus there exist low-energy observables in which the quantum effects 
induced by virtual heavy particles do not vanish when the mass of these particles goes to 
infinity. 

Most prominent among the non-decoupling effects is the p parameter [2], which pro- 
vides a measure of the relative strength of neutral and charged current interactions in four 
fermion processes at zero momentum transfer. At tree level p = 1 due to a global acciden- 
tal SU{2) symmetry, the so called custodial symmetry, p can receive radiative corrections 
only by those sectors of the SM which break explicitly the custodial symmetry, namely the 
hypercharge and the Yukawa couplings which give a different mass to the components of 
fermion doublets. In this latter case, the contribution to the p parameter is proportional to 
the mass splitting, therefore the leading contribution comes from the top-bottom doublet. 

At one loop, the p parameter has a quadratic dependence on the top quark mass, 
A/^ 1 ) ~ Gi?7?ij, and a logarithmic dependence on the Higgs mass, Ap*- 1 ) ~ g' 2 log 
Two loop corrections at the leading order, i.e. Ap^ ~ G 2 F m\, and at the next-to-leading 
order, i.e. A// 2 ) ~ G 2 F m 2 M§, in the top quark mass were computed in the limits mn — > 
and ran 3> m t in Refs.[3l 0J and for arbitrary Higgs mass in Ref.jS]. It turned out that 
due to accidental cancellations, the subleading corrections at two loop are larger than the 
leading ones [6]. At three loops, the computation of the leading top quark corrections, 
Ap^ « Gpmf, in the massless Higgs limit, was carried out in Ref.[7]. The complete 
dependence on the Higgs mass at three loops was obtained in Ref. [8]. Numerically it 
was found that this contribution to Ap is quite large and provides a sizeable correction 
(~ 36%) to the leading electroweak correction at two loops. However, the size of the 
three loop correction is only about 2% of the much larger two-loop subleading electroweak 
correction. Moreover, the perturbative series of the leading top quark contributions to the 
p parameter has alternating signs up to three loops. 

This raises the issue of the convergence of the perturbative expansion (it might be that 
this series is divergent, but Borel summable) and calls for a better understanding of higher 
order radiative corrections. It would be highly desirable to have a simplified framework 
in which the leading top quark contributions to the p parameter can be computed to all 
orders in perturbation theory and eventually summed up. The actual calculation of the 
leading radiative corrections in the top quark mass is greatly simplified by the observation 
that to obtain them it is enough to consider the lagrangian of the SM in the limit of 
vanishing gauge coupling constants g,g' — > [3J. This gaugeless limit provides an efficient 
way of reducing the number of Feynman diagrams to be computed and it has been used 
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in the two and three loop computations mentioned above. 

The SU(Nf) x [7(1) electroweak model in the large iVp-limit [HI [TO] provides an ideal 
framework in order to compute the leading top quark contributions to the p parameter 
to all orders in perturbation theory. In fact, since only one-loop graphs contribute to 
the top quark self-energy at the leading order in the large iV^-limit, the exact top quark 
propagator can be obtained simply by resumming one-loop self-energy insertions. In this 
way, one takes into account the finite width effects due to the fact that the top quark is an 
unstable particle. The resulting Dyson propagator contains, in addition to the complex 
pole corresponding to the top quark, a tachyon pole in the euclidean region, p = -A 2 T , 
which spoils causality and makes all the Wick-rotated Feynman integrals ill-defined. 

Since the tachyon pole is a non-perturbative effect one can still compute the p parame- 
ter by using the resummed top propagator instead of the Born one, provided an expansion 
in powers of the top Yukawa coupling is taken before performing the Wick rotation. All 
the coefficients of this perturbative expansion can be computed analytically. It turns out 
that they are all positive and grow factorially with the number of loops, thus their power 
series is divergent and not Borel summable. 

In order to go beyond the perturbative approximation, one has to regularize the inte- 
grals containing the resummed top quark propagator. In this connection, the introduction 
of an UV cutoff at A < At has been proposed in Ref.[TT]. However, this procedure breaks 
gauge invariance. We have adopted a different strategy devised in Ref. [15]. Assuming 
that the occurrence of the tachyon pole is not due to the inconsistency of the theory 
under consideration, but of the intermediary expansion technique used, it is reasonable 
to circumvent the ill-defined part by an adequate subtraction of the tachyonic pole. In 
particular, we have chosen to subtract the tachyon minimally at its pole from the exact 
top propagator. One should be careful in doing this because the tachyon pole contributes 
to the Kallen-Lehmann spectral function. The correct normalization of the latter, which is 
crucial in order to guarantee renormalizability, is recovered after the tachyonic subtraction 
by a suitable rescaling of the top propagator. 

This procedure allows us to find a tachyon-free representation of the exact leading top 
contribution to the p parameter which can be estimated numerically and compared with 
the expansion of Ap at any fixed order in perturbation theory. 

The paper is organized as follows. In Sect. [5Jthe SU(Nf) x U(l) model in the large 
./V^-limit is presented and its spectrum is briefly discussed. The top quark self-energy at 
the leading order in Np is computed in Sect. The subtraction term for the removal 
of the tachyonic pole is discussed in Sect. 13.11 In Sect. 0] we check the validity of Ward 
identities connecting the self-energies of gauge bosons and of unphysical scalar particles 
computed with the resummed top propagator. The leading top quark contribution to the 
p parameter is computed to all orders in perturbation theory in Sect. [5] In Sect. [6] the 
tachyon-free representation of the resummed top propagator is used in order to compute 
nonperturbatively the leading top contribution to the p parameter. Finally the conclusions 
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are given in Sect. [7J In the Appendices we have collected technical details concerning the 
computation of the radiative corrections to the p parameter. 



2 The model 

In this Section we shall consider a SU (Np) x (7(1) gauge theory in the large iVj?-limit. For 
the sake of simplicity only, one generation of fermions will be taken into account. The non- 
abelian gauge fields W® transform according to the adjoint representation of SU (Np), thus 
a = 1,2 ... , (Np — 1). The abelian gauge field is denoted by B^. Left handed fermions 
transform according to the fundamental representation of SU(Np), therefore they are 
sorted in iVp-plets, 

Ll -"- L -( ) • « L - = ( ) ■ i = »■ 2 - ■ ■ ■ N ? ~ 1 • « 

while right handed fermions are singlets as in the SM 

if = oj+lj , t R = uj+t , bf = uj + bj , j = 1, 2, . . . N F - 1 , (2) 

where co± = 1= ^ 75 are the chiral projectors. All fermions are taken to be massless except 
for the top quark. Finally, the scalar sector of the model is a gauged 0(2Np) linear sigma 
model in the broken phase. The scalar fields are sorted in a complex iV^-plet 

•=(^ + ffi > + ix(.»)' J ' = 1 ^- W '- 1 ' (3) 

which transforms according to the fundamental representation of SU(Np). The classical 
lagrangian of the model is given by 

c = -- A (d^wi - d v w; + gf ahc K w i) 2 ~ 7 ( d » B - - d » B » 



+(D^) ] (D^^) - 



+L L iJp L L + Q L iJp Q L + lj R iJp if + t R ip t R + bf iJp bf 

+ -L yt (v + H)it - ^=y tX tj 5 t - y t (f)J bf t L - y t 0+ t L b] , (4) 

where f abc are the structure constants of SU(Np), g and g 1 are the SU(Np) and U(l) 
gauge coupling constant respectively, yt is the Yukawa coupling constant of the top quark, 
A and v are the quartic coupling constant and the v.e.v. of the Higgs A^-plet respectively, 
while is the covariant derivative for matter fields. 

The spectrum of the theory has been discussed in Ref. [9j [10]. For our purposes 
it is enough to mention here that there are Np — 1 copies of the W boson and of the 
charged unphysical scalar field, <f>. In the quark iVp-plet the upper component, which we 
can identify with the top quark, is massive, while the remaining Np — 1 components are 
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Figure 1: Enhanced one- loop self-energy graphs for the top quark. 

massless. These latter are essentially copies of the bottom quark. The massless lepton 
Np-plet does not play a role in the following. 

In order to perform the large N F limit we make the following assumptions: 

g 2 N F , g ,2 N F , y 2 t N F , XN F , ^- (5) 

will be held fixed while taking N F to infinity. This amounts to say that the mass of every 
particle in the SM is of order O(l) in the large iV^-limit. For the sake of simplicity in 
what follows we shall consider a SU(1 + N F ) x U(l) gauge theory. In this way the SM 
corresponds to the case N F = 1. 



3 Top quark self-energy 

In this Section we shall discuss the one-loop top quark self-energy in the SU(1 + N F ) x U (1) 
model and its renormalization in the on shell scheme. Moreover, the subtraction of the 
tachyonic pole from the renormalized top propagator will be presented. 

Due to the presence of N F copies of the W boson, the unphysical charged scalar 4> and 
the b quark in the model under consideration, two graphs contributing to the top quark 
self-energy, St(p), are enhanced by a factor N F . These graphs are depicted in Fig.([T]). The 
remaining one- loop self-energy graphs are of order 0(1). Moreover it is straightforward 
to prove that higher order 1-PI graphs give subleading contribution to the top quark 
self-energy in the large A^r-limit. This means that the computation of the two graphs in 
Fig.([T]) is enough to get the exact top quark self-energy at the leading order in the large 
iVp-limit. 

By direct inspection it turns out that graph (1) is proportional to the squared mass of 
the top quark, while graph (2) is proportional to the squared mass of the W gauge boson. 
We are interested in the case in which the top quark is the heaviest (and only) mass scale, 
thus we neglect the contribution of graph (2). The dimensionally regularized contribution 
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of graph (1) to £t(p) is given by 

^t(p) = G F m 2 tfi N F -^-j + log (-j2~ »e) i>^+ , (6) 

B 

where m tj o = -^ytv is the bare top quark mass, while A# is a regulator dependent 
quantity with the dimension of a mass whose explicit expression is not needed in the 
subsequent analysis. 

The natural way to incorporate the finite width effects of an unstable particle is the 
resummation of the corresponding self-energy insertions. This leads us to consider the 
Dyson resummed top quark propagator instead of the Born one. We report here only the 
component of the resummed propagator with positive chirality because it is the only one 
which gives a non vanishing contribution to the p parameter. 

D t {p) = 2— ' where 

ao[P )P ~ mfo 



V2 



a Q {p 2 ) = l-—-,G F m 2 N F ^_ + log (- ^ - * e ) . (7) 



16vr 2 



2 / P 



D - 4 ' 6 V A| 



2 



It should be noted that the above expression is valid to all orders in the interaction 
strength and provides the exact bare top quark propagator at the leading order in the 
large iV^-limit. 

In order to renormalize the top quark self-energy, we require that the real part of the 
denominator of the Dyson propagator in Eq.([7|) vanishes when computed at p 2 = m 2 . This 
allows us to express the bare top quark mass in terms of the subtraction point mt 

r 2 / Tn 2 \ i 

m t,o = Re[a (m t 2 )] m 2 t = m 2 - a t [jj^ + log (j^J \ m t , (8) 

where we have introduced a shorthand notation: at = j0? G F m 2 N F . 

By substituting the above equation into Eq.([7]), we obtain the on shell renormalized 
top quark propagator at the leading order in the large- N F limit 

D t {p) = , i^^, . , where a{p 2 ) = 1 - a t log (- 4 - ■ (9) 



a \P)P — + ie V mi 

Notice that in order to achieve the above result a suitable wave function renormalization 
for the left- and right-handed components of the top quark field must be imposed. On the 
top quark mass shell, p 2 ~ m 2 , we have 



A(p) = — 2~- F" ' ^ 10 ) 

p z — mi + imtit 



where Tj = tt at mt is the total decay width of the top quark. We remark that this is 
the exact width of the top quark at the leading order in the large A^ir-limit in the narrow 
width approximation, i.e. for Tt <C mt- 
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3.1 Tachyonic regularization 

Besides the complex pole corresponding to the unstable top quark, the propagator in 
Eq.(|9|) contains a tachyon pole. Its euclidean position, p 2 = — A 2 ^, can be obtained by 
solving numerically the following equation 

1 A 2 
1 + -y = a t log (A T ) , where \? = -\ . (11) 

Xrp Ulj. 

The tachyon pole induces causality violation effects in the theory and makes all the Wick- 
rotated Feynman integrals ill-defined. Thus in order to obtain sensible results, it must be 
removed from the integrals involving the top propagator. 

One simple and consistent way to deal with the tachyon is to regard the model under 
consideration as a low-energy effective theory and to introduce explicitly a cutoff under 
the tachyon scale, A < At, [11] ■ However, this procedure breaks gauge invariance and 
thus we prefer to use a different approach. Assuming that the tachyon is a mere artifact 
of perturbation theory, we modify the top propagator in Eq.([9]) by subtracting minimally 
from it the tachyonic pole. In order to determine the correct normalization of this tachyon- 
free representation of the top propagator, let us consider the spectral representation of the 
resummed top propagator @ 

D t (p) = ip'uj + / ds — , where p(s) = p T (s) + p, (s) 0(s) . (12) 

J-oo P ~ s + ie 

Notice that due to the tachyonic contribution to the spectral function, 

p T (s) = k5(s + A T ) , where k= — -5- ~ — exp ( — — V (13) 

1 + a t A T a t V at J 

is the residuum at the tachyon pole, the exact top propagator Q does not satisfy the 
usual Kallen-Lehmann representation. The other contribution to the spectral function, 
which comes from the positive part of the spectrum, is given by 



p+(s) = ^ -\ -2 r • 

a t A log (A) -1 



1 mi 



H 

Clearly the removal of the tachyonic pole is necessary in order to find an expression 
for the resummed top quark propagator which respects causality and satisfies the Kallen- 
Lehmann representation. On the other hand, the contribution of the tachyon pole is 
crucial in order to ensure the correct normalization of the spectral function, as can be 
easily checked numerically 

/+00 r+00 
ds p(s) = k + / ds p + (s) = 1 . (15) 
-00 Jo 

In order to compensate the tachyonic contribution to the spectral function a suitable 
rescaling of the top quark propagator should be performed. This leads us to the following 
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tachyon-free representation of the exact top quark propagator 

n m - i ^ >UJ+ \ a ^ — 1 (Mi\ 

t[P) 1 - K U(p 2 ) p 2 - ml p 2 + A 2 T \' { ' 

Some comments are in order, i) Since the subtraction term vanishes to all orders in 
perturbation theory, this tachyonic regularization can be regarded as a prescription for 
summing up the perturbative expansion in which a tachyon does not show up in the 
theory's spectrum, ii) Although not unique, our prescription preserves gauge invariance 
to all orders in perturbation theory, hi) The prefactor 1/(1 — k), which ensures the correct 
normalization of the spectral function after the subtraction of the tachyon pole, is crucial 
in order to prevent the appearance of spurious divergences in the radiative corrections to 
the p parameter. 

4 Ward identities 

In this Section we present the exact top quark contribution to the one-loop vector (i.e. Z 
and W) and scalar (i.e. x an d 4>) self-energies at the leading order in the large A^-limit. 
Since, in general, resummation of self-energy insertions can spoil gauge invariance, we will 
also check the validity of the Ward identities connecting these self-energies. 

It is convenient to perform the computation of vector and scalar self-energies in a pure 
SU(1 + Np) gauge theory without hypercharge. This simplifies the actual computation 
without modifying the leading top quark contribution to the p parameter. 

Let us recall the expression of the sine and cosine of the Weinberg angle and of the 
electric charge in the zero hypercharge limit, g' = 0. 

v g 2 + g' 2 Vg 2 + g 2 Vg 2 + g 2 

As a consequence of the above relations, as expected, M\y = Mz- 

The coupling of the Z boson to the top quark in the SM is given by 

£1^7^- +5 t + w+). (18) 

If we switch off the hypercharge, we find 

f 2 s w 1 1 \ g + 2 es w 

eg t =e[-- + - =-, eg+ = -- = 0. (19) 

V 6 Cw 2 SwCw' 2 6 Cw 

The couplings of the W boson to fermions do not change if we set g' = 0. 

For our purposes it will be enough to consider only those graphs contributing to the 
vector and scalar self-energies where at least one virtual top quark is exchanged. At the 
one-loop level this task amounts to the computation of two Feynman graphs which are 
shown in Fig. ([2]) for the case of the vector self-energies (the other graphs are the same 
with scalar external legs). The required amplitudes can be easily obtained by using the 
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Figure 2: Top contribution to the Z and W boson self-energies at one loop. The solid line 
with a bubble denotes the resummed top propagator. 



SM Feynman rules for vertices in the zero hypercharge limit (see Eq. (|19|) ). the resummed 
top propagator in Eq.© and the Born propagator of a massless b quark. In order to avoid 
cumbersome expressions and since they do not play any role in the proof of the validity of 
Ward identities, we do not show in this Section the contributions coming from the tachyon 
subtraction term. Their discussion is deferred to Sect. [6j 
The Z self-energy in the zero hypercharge limit reads 



Ef(p 2 ) 



■g 2 N c 



d°q a(q 2 )a((q - p) 2 )[g^{q 2 - q • p) - 2q»q v + q^p u + q v p»\ 



(2vr 



[a{q 2 )q 2 - m 2 ][a((q - p) 2 )(q - p) - m\ 



(20) 



where N c is the number of colours. 
The x self-energy, whose expression does not depend on the hypercharge g' , is given by 



W Nc ~ 

M 2 
1V1 W 



a((q-p) 2 ) Q-{Q-p)-m 2 
(2ir) D [a(q 2 )q 2 — mf ] [a((q — p) 2 )(q — p) 2 — m 2 ] 



d u q 



(21) 



By dotting p^ , p u in Eq. (f2"Uj) we get 

^tJl\ _ i 9 2 N c f d D q a{q 2 )a{{q - p) 2 )[q 2 p 2 + {{q - p) 2 - q 2 ) q ■ p] 



PfiPu^z'yP 



(27T) D [ (,2) q 2 _ m 2] _ p )2)( ? _ p y _ m 2] 



i g 2 N c m 2 



a{{q-p?) {p 2 



P) 



(2*) D [a(q 2 )q 2 -m 2 ][a(( q -p) 2 )(q- P y 



(22) 



Notice that in Eq. (|2ip there is a term proportional to mf which is absent in the second 
line of Eq. (j22p . This term, however, cancels if one takes into account the contribution of 
Higgs tadpoles to the scalar self-energy S x . 



->tad 



-ig 2 N t 



M 2 W 



1 



d D q 

(2vr) D a(q 2 ) q 2 - m 2 



(23) 



If we add S* ad to the self-energy in Eq. (j2ip we get 



S x (p 2 ) = ig 2 N c 



ml 



d D q 



a((q ~ p) 2 ) (q ~ p) -p 



M^J (27T)D [aiq 2 )q 2_ m 2 ][a{{q _ p) 2 ){q _ p) 2_ rn 2 ] 



(24) 
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By comparing Eqs. (|22|) and (|24|) we find 

^^V) + £ x (p 2 ) = 0. (25) 

We consider now the one-loop W self-energy. 
The one-loop <f> self-energy is given by 

V I "2\ ■ 2 at m t f dD <l Q-(q-p) , ov x 

By dotting p^ , p u in Eq. (f26j) we get 

d D q a{q 2 ) q 2 p 2 + a(q 2 ) [(q - p) 2 - q 2 ]q • p 



(2vr) D [a(q 2 )q 2 -m 2 }(q-pY 

i f d D q p 2 -q-p 

(2ir) D [ ( g 2) q 2 _ m 2] ( ? _ p )5 



% 9 N c™ t / 7 ^- w — — — — (28) 



By comparing Eqs. (|28p and ()27[) one finds that a relation analogous to the one in Eq. (|25|) 
does not hold. However if one adds the contribution of Higgs tadpoles 

Er ._ i/W ^/^_L_, (29) 

to the self-energy in Eq. (|2"T]) one gets 

E^(p ) = i 5 iV c —2- / — — n — — — — -2 . (30) 

M W J ( 27r ) Hq 2 ) q 2 - mf\{q - p) 

Now it is straightforward to prove the following Ward identity 

^^(P 2 ) + S,(p 2 ) = 0. (31) 
1Vi w 

5 Perturbative top contributions to the p parameter 

In this Section we shall use the resummed top propagator in Eq.© in order to compute the 
leading contributions in the top quark mass to the p parameter. It turns out that due to the 
presence of a tachyonic pole, the resulting expression for Ap is ill-defined. However, while 
the tachyon is a nonperturbative effect, all the coefficients of the perturbative expansion 
of Ap in powers of the interaction strength at are well-defined and can be computed 
analytically. 

The p parameter is usually defined as the ratio between the neutral and charged current 
coupling constants at zero momentum transfer 

P ~ Jcc(0)- 1-Ap- (32) 
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-^Cc(O) is given by the Fermi coupling constant, Gf, determined from the //-decay rate, 
while Jnc(Q) can be measured in neutrino scattering on electrons or hadrons. Notice that 
this definition of the p parameter is process dependent since, in general, the radiative 
corrections depend on the hypercharge of the particles involved in the scattering process. 
However, the leading contributions in the top quark mass to Ap are universal. 

M 2 

At tree-level the p parameter is given by p = . f - 2 w :i =1. At the leading order 

M z c w 

in the top quark mass radiative corrections to p stem from the transversal parts of the 
(unrenormalized) self-energies of the vector bosons W and Z 

Ap = — T o- , wnere 

^(p 2 = 0) = U v g^,V = W,Z . (33) 
Hz and ILw can be obtained by setting p = in Eqs. ([20|) . ([26]) respectively. 
TT _ \2 N U 2 X fd D q a 2 (g 2 )g 2 

^-^N c (l-l)j^ [a{q2)q :l^ ]q2 . (35) 

By substituting Eqs. (j34]) and (|35l) into Eq. (j33l) . it is now straightforward to derive the 
leading top quark contribution to the p parameter. Since it turns out that the result is 
both IR- and UV-convergent we can compute it directly in D = 4 % 

AP = I ^ Nc M% J [atf) - m?] 2 g * • (36) 

In the free field case {a% = 0), the above equation gives the well known result 

g 2 N c m? N C V2 2 N c 
= ~64tt 2 ~ = ~IQtt 2 ~ m * = iV> ' (37) 

After expanding the denominator in Eq. (|36p about at = 0, one can perform a Wick 
rotation and integrate over the solid angle 

n f ~ Jo (i + x) u+2) 



j=0 



where we have introduced the dimensionless variable x = — ■ 
Let us consider a slightly more general class of integrals 



Ii(j) = (j + 1) dx- ^,withi,leN,i>I + l. (39) 

Jo (l + x) u + ~> 



3 Notice, however, that the integral in Eq. (|36l) is ill-defined due to the tachyonic pole present in the 
resummed top propagator @. As a consequence, the Wick rotation cannot be performed because the 
resulting integral would be divergent. 
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Clearly, we are interested in Io(j), but if we integrate by parts the latter, we immediately 
find a linear combination of Iq(J — 1) and (see the identity in Eq. (|40p ). The case 

j = I must be dealt with separately because boundary terms contribute. Integration by 
parts can be applied to giving rise to the following recurrence relation 

m = f±M w . l) + U±M^ lM(j) . (40) 

Eq. (|4U|) allows us to express Io(j) in closed form by means of a linear combination of 
simpler integrals 

3 poo / i \^ 

I O (i)=£ c <(iH (2) (0, where l\ 2 \l) = / dx . (41) 

i =0 J o (1 + x) 

The coefficients c/(j) are related to the combinatorial problem of grouping together I 
objects out of a total of j without repetions. Their explicit expression and some useful 
properties are reported in Appendix [Al 

(2) 

In order to compute I, (/), let us consider the following change of variable: logx = t 

/+oo 1 1 

Since the integrand in Eq. (|42|) is an odd function if I is odd and an even function otherwise 
we immediately obtain the following result 

4? (20 = J dt et + e _ t + 2 , ISU21 + 1) = . (43) 



The integrals in Eq.(j43J) can be computed analytically by exploiting the properties of 
polylogarithms, as shown in Appendix [Bj We report here only the final result 

4 2) (2Z) = 2 (2/)! (l - 2^- 2l) ) ((21) . (44) 



Putting all the pieces together we obtain 

00 [j/2 

N F 



A P=§ L Y, r i a ¥ +1) , where r, = ^2c 2/ (i)(2/)!(l-2( 1 - 2i ))c(2/), (45) 
7 i=o 1=0 



where with [ • ] we denote the integer part of a real number. As an example, we report 
here the perturbative expansion of the leading top contribution to the p parameter up to 
terms of order five in at 



Nc 
N F 



f ^ 2 )a? + (l + ^) a 3 + ( 1 + ^ 7r 2 + ^ 7r 



(46) 



By making use of the asymptotic estimates of the combinatorial coefficients in Eqs. dTOD 
and (|73p we can easily find the leading order behaviour of the coefficients of the perturba- 
tive expansion of the p parameter. In particular, it turns out that rj ~ (j + 1)! for j ^> 1. 
Thus the perturbative expansion of the p parameter is factorially divergent and not Borel 
summable, being a fixed sign power series. 
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5.1 Wave function renormalization of unphysical scalars 

In order to perform a cross check of our result in Eq. (|45p . we compute the radiative 
corrections to the p parameter by means of the wave function renormalization of the 
unphysical scalar fields. This latter approach is related to the one pursued in Sect. [5] by 
gauge invar iance. 

Let us consider the kinetic terms of the scalar part of the SM lagrangian. The UV 
divergences that show up in radiative corrections can be reabsorbed by introducing suitable 
wave function renormalization constants in the following way 



2 



+ 



gv 



Z u ) + other terms . 



2 V p/v 2c w 

The renormalized masses of the gauge bosons are given by M\y = \[Z§ %f and Mz 



(47) 



since Zg 



1 + ^ 



, S = d 



p 2 =0 



p 2 =0 

X ■ 



(48) 



The self-energies of the unphysical scalars have been computed in Sect. 0] (see Eqs.([2 
([30]) ). The p parameter reads 

A a r d D q 



Ap = iq 2 N c — I s /- 



q-p 



(49) 



(2tt) d [ a (?) f - m 2 ][a((q - pf) (q - pf - mf] q 2 ' 

We now develop the denominator about p = 0. For the computation of the derivative 
w.r. t. p 2 it is enough to keep terms of the order of q ■ p, because a term proportional to 
p^ is already present in the numerator 



9 » T mf d 



d D q (q-p)\a{q 2 )-a t ) 



(27r) [ a tf 



m 



2l3o 



(50) 



The tensor reduction can be performed immediately. In fact, since the denominator does 
not depend on p, q^q u must be proportional to g^ u . Moreover we can work in four 
dimensions because the above integral is both IR- and UV-convergent. 



i 2 v4 fA jWU 

2 y C M&J (2tt) 4 [ a (qi) q 2 



■9 N c 
2 y M? 



mi 



m 



a t q* 



2l 3 r,2 



d A q 

(2^) 4 [a(q*) q* - m%] 3 , 



+ 



d 4 q : 

[a(q 2 ) q 2 



m 



• (51) 



Notice that the second term in the last line of the above equation is simply 2Ap, as can 
be seen from Eq. (|36p . thus finally we are left with 



Ap = - g 2 N c 



ml 



d 4 q 



mr 



m 



213^2 



(52) 



4 Notice that this works because the derivative w.r.t. p 2 is both IR- and UV-convergent, otherwise one 
looses finite parts by computing the derivative in this way. 
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After expanding about at = 0, one can Wick rotate the above expression and integrate 
over the solid angle, finding 



N c v^, w , i f°° (1 + a 
Ap=^«* J>- + 2)(j + lK / dx K - 

F j =Q JO (1 



1 + at x)x> (logx) J 



(1 + x) 



(3+3) 



(53) 



At the leading order in the interaction strength ctj, we immediately find the result in 
Eq.([37D 

2 N r . 



A N * 



/ dx — 
Jo (1 



(i + x y n f 



a t . 



(54) 



Let us consider the coefficient of a\ , with j > 1 (for the sake of brevity in the following 
equations the common prefactor M^- at is omitted) 



dx p]ogx^_ + (j + i)j 



(1 + x) 



00 , x^-^riogx)^ 

ax 



(1 + x) 



(i+i) 



-(j + i)i / dx 



00 ^0- 1 )(logx) (jf ~ 1 ' ) 



(55) 



(l + xf +2) ' 

Notice that the second term in the first line of the above equation is given by (j+1) Io(j— 1). 
After integrating by parts the first term in the first line of Eq. (|55p . we find 

(3 + 1)3 dx- L^- + ( J + l)/ o(j -l). (56) 

Jo (l + x) u+) 

By applying iteratively integration by parts on the first term of the above equation and 
by using Eq. (j4"T|) , we get 



l>0-i)(j-l) 



1=1 



I + 1)1^(1-1] 



(57) 



(3) 

where the integrals L (I) are given by 



7 f> (()=2 r dI (i!54 =2 rv 

Jo (1 + x) J -00 



e* t l 



(1+*) 



3 • 



(58) 



(1 + x)- 

In the above equation a change of variable t = log x has been performed. We integrate by 
parts the integrals in Eq. (|58p and we postpone the evaluation of the boundary contribu- 
tions (for the notations see Appendix [B]) 



2 dt 



e* t l 



(1 + e*) 3 " (1 + e*) 2 



+ t-l dt 



e t t (i-i) 



(l + e*r 



I dt- 



1 + e* 



t i [Li_ 1 (-e < ) + Li (-e*)]+/r dtt^U-i(- e*)+ dt Li ( - e*)l . (59) 



Some comments are in order, i) The term — t'LL_i( — e*) vanishes when evaluated at the 
boundaries, i.e. ±00, therefore it can be neglected, ii) The first integral in the second line 
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of the above equation vanishes when I is even, while for I odd its result is given in Eq. (|79|) . 
iii) The last integral in the second line of the above equation is divergent if evaluated at 
(positive) infinity, however this divergence is exactly compensated by — i Lio( — e l ). We 
compute this latter integral by exploiting the properties of the poly logarithms. 

I r-l 



[dtt^Uoi-e*) = X>) r+1 [fl(l-sj\ i^>Li r (-e*) 

J r=l s=0 

= - Jdtt 1 Li_i( - e*) + t l Li ( - e*) . 



(60) 



The second line of the above equation follows immediately by a comparison between the 
sum in the first line of the same equation and Eq. (|79p . Putting all the pieces together, we 
find 



7f )(Z) = - 1) + lj; 2 \l) for I > 1. 



If we plug the above result into Eq. (|57j) . we get 



j 



In the above equation we have used the identity in Eq. ([68l) . 



(61) 



(62) 



6 Nonperturbative top contribution to the p parameter 



In this Section we shall use the tachyon-free representation of the resummed top propagator 
(|16l) in order to compute nonperturbatively the exact leading top quark contribution to 
the p parameter at the leading order in the large A^-limit. 

The contribution of the tachyonic subtraction term in Eq. (|16p to the one- loop self- 
energies of the W and Z vector bosons at zero external momentum can be easily computed. 
We show here the results. 



II; 



9 Z N C (1 



-)- 

DJ (1 



1 



d D 



q 



a 2 (Q 2 )- 



(2tt) d 1 [ a ( g 2 



m 



21 2 



+ 



2k 



k 2 A 2 t 



2k m 2 



( g 2 + A 2 )2 [a{q 2 )q 2 



m 2 ](q 2 +A 2 T ) 



(63) 



n 



w 



ig 2 N c (i-^-) 



d D C 



mi 



D/ 1 — k J (2tt) d I [a(q 2 )q 2 



+ 



kA 2 -, 



q 2 (q 2 + A 2 T ) q 2 



(64) 



By using Eqs. (j63j) and (JM|) into Eq. (j33j) and by setting D = 4, one can write down a 
tachyon-free representation of the leading top contribution to the p parameter 

2 



Ap 



4 y M? 



1 



(1 



d"q 1 



(2vr) 4 q 2 



1 



a(q 2 ) q 2 



+ 



mi 



K X 2 ^ 

2 + 



(65) 



15 



Notice that the perturbative expansion of the above result coincide order by order in 
at with the factorially divergent and not Borel summable perturbative series in Eq. (l38D 
since the additional terms, proportional to k, vanish to all orders in perturbation theory. 
However, the integral in Eq. (|65p is now well-defined and thus the Wick-rotation and the 
integration over the solid angle can be performed directly on it. The result of these 

2 

operations can be expressed in terms of the dimensionless variable x = — ^ as follows 

2 



Ap 



— — at — 

N F *(i 



1 



dx 



1 



a{— x) x + 1 



K Xj, 

v — A^ 



(66) 



The above integral can be computed numerically (for instance with the help of Mathemat- 
ica) for an arbitrary value of the interaction strength at, allowing us to make a comparison 
between the exact nonperturbative result and its perturbative approximation at any fixed 
order in at- 

In Table [T] we show the position of the tachyonic pole (divided by mt), the residuum 
at the tachyon pole, n, and the leading top contribution to the p parameter (omitting the 
prefactor N c /Np) for some values of at- 



a t 


Ay 


K 


Ap 


0.02 


7.2 • 10 10 


9.6- 10~ 21 


0.021 


0.04 


2.7- 10 5 


3.5 • 10" 10 


0.042 


0.06 


4.2 • 10 3 


9.6 • 10" 7 


0.065 


0.08 


518.0 


4.7 • 10" 5 


0.090 


0.10 


148.4 


4.5 • 10" 4 


0.118 


0.20 


12.38 


0.032 


0.249 


0.40 


3.805 


0.147 


0.329 


0.60 


2.602 


0.198 


0.341 


0.80 


2.141 


0.214 


0.344 


1.00 


1.895 


0.218 


0.344 



Table 1: Numerical values for the tachyonic pole, its residuum and the leading top contri- 
bution to the p parameter. 



The exact numerical result, Ap(at), shows a typical saturation behaviour (see Fig. [3]) 
for at > 0.2 which cannot be reproduced by the perturbative expansion of the p param- 
eter (|45p at any fixed order since all the expansion coefficients are positive (see Fig. U]). 
However for small enough values of the interaction strength, say at < 0.2, the agreement 
between the nonperturbative exact result and its perturbative approximation (starting 
with terms of order 0(af)) is very good. Finally, it should be noted that since the per- 
turbative expansion of the p parameter is a divergent asymptotic series, the perturbative 
approximation of the exact result can be improved by adding further terms to the series 
only up to a certain order, beyond which the approximation gets worse and worse. 



16 




Figure 3: Exact leading top contribution to the p parameter 
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7 Conclusions 



In this paper the SU(Np) x U(l) model at the leading order in the large A^-limit has been 
used in order to compute the exact leading top quark contribution to the p parameter and 
its perturbative expansion to all orders in the interaction strength at- 

Since only one-loop graphs contribute to the top quark self-energy at the leading 
order in the large A^-limit, the exact top quark propagator can be obtained simply by 
resumming one-loop self-energy insertions. In this way, one takes into account the finite 
width effects due to the fact that the top quark is an unstable particle. On the other hand, 
this Dyson resummed propagator contains a tachyon pole in the euclidean region which 
spoils causality and makes all the Wick-rotated integrals ill-defined. We have regularized 
the resummed propagator by subtracting the tachyon minimally at its pole. Although 
not unique, this procedure allows to define a tachyon-free representation of the exact top 
propagator which respects gauge invariance. 

The validity of the Ward identities connecting the self-energies of vector bosons and 
of unphysical scalar particles, computed by using the resummed top propagator instead 
of the Born one, have been checked. These vector and scalar self-energies then have been 
used in order to compute the leading top contribution to the p parameter in two different 
ways as a further check of gauge invariance. It turns out that the perturbative expansion 
in powers of the interaction strength at of the p parameter is factorially divergent and not 
Borel summable. 

However, after having subtracted consistently the tachyonic pole the expression for the 
leading top contribution to the p parameter can be evaluated numerically and compared 
with its perturbative approximation. The agreement between the exact result and its 
perturbative expansion (starting with terms of order 0(af)) is very good for at < 0.2 which 
in the SM, i.e. for Np = 1, corresponds to a top quark mass of 1.4 TeV. Moreover, the 
exact numerical result shows a typical saturation behaviour which cannot be reproduced 
by the perturbative expansion of the p parameter at any fixed order, since all the expansion 
coefficients are positive. 

Though the subtraction of the tachyon pole is determined by the demand of causality, 
the procedure is not quite unique, since the correction factor that is needed to insure a 
properly normalized spectral density could have been different from a constant. However, 
one can consider this correction factor, which is given by the residuum of the tachyon 
pole, as an estimate for the uncertainty in the calculation due to non-perturbative effects 
or effects of new physics at high energy. The uncertainty is at most of the order of 20%. 
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A Combinatorial coefficients 



In this Appendix we give the explicit expression of the combinatorial coefficients q(j') 
introduced in Eq. (|4ip and we show some of their properties. 

The recurrence relation in Eq. ljiT)]) can be applied to the integrals h{j) (see Eq. (j39j) ) as 
long as I < j. Thus, starting from Io(j) and applying repeteadly the recurrence relation, 
one ends up with a linear combination of /;(/), with < I < j. The coefficients of this 
linear combination are 

co(j) = 1 , 

1 l j-l+i 

Q(j) = ^n E r i5 with n = l,2,...i-Z + l. (67) 

i=lr,=r(i-i)+l 

By using the definition of ci(j) in Eq. (l67p . it is straightforward to show that Cj(j) = 1. 
Moreover, another useful relation which can be easily proven is the following 

q(i)-Q(i-l) = yC (z _ 1) (j-l). (68) 

The knowledge of the asymptotic behaviour of q(j) for j>l will allow us to determine 
the large order behaviour of the perturbative expansion of the p parameter in Eq. ()45D , 
For this purpose, by making use of the following identity 

Z*='£*-'ZV*'—-— , (69) 

k=r+l k=l k=l 

we eventually find 

ci(j) = \ (y|) 2 + O0' (,_1) ) , for j » 1 , / « J . (70) 

In order to obtain an asymptotic estimate of q(j) which holds for / ~ j, it is convenient 
to write down an expression for the 'last' coefficients C(j_n(j) 

*y-oc*)=[rio--*+i)]n J £ v- (n) 

k=l i=l rj=r (i _ 1) +l 1 

The sums of the reciprocals of natural numbers in the above equation can be rewritten in 
terms of products of the harmonic numbers 

/! n e \. - ( h ^ 1 > where = e \ ■ ( ?2 ) 

i=lri=r^^ r) +l 1 r=l 

Since H(j) ~ logj for j 3> 1, the leading order asymptotic behaviour of C(j_^(j) is given 
by 

c i-^') = j] ( J L i)i ( lo §^' + °(( , for j » 1 , 1 « j . (73) 
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B Integrals 



In this Appendix we compute the integrals in Eq. (|43p by making use of the properties of 
poly logarithmic functions. 

The polylogarithm Li s (z) is, in general, a special function defined by the following 
series 

OO f. 

Li a (2) = ^2 ' Vz ' s G C ' with \ Z \ <1 - ( 74 ) 

k=l 

By analytic continuation it is possible to extend the domain of the polylogarithm over a 
larger range of z. Notice that for some values of the parameter s, it is possible to express 
the polylogarithm by using elementary functions. For instance 
00 00 
U Q (z) = J2* k = 1 J t- z , Lh(z) = ^2 Z T = -log(l-z). (75) 
k=i fc=i 

By using the definition in Eq. (|74p and by integrating the series term by term it is straight- 
forward to prove that 

U s+1 (z) = fdt ^ , thus j- U s+1 (z) = ^ . (76) 
Jo t dz z 

We list here some properties of the polylogarithms which are needed for the computation 
of the above mentioned integrals. 

lim Lis iz) = . 

\z\-yQ 



Li 



:(-!) = - (l - 2 [1 ~ s) ) Lis(l) = -(l - 2^)) C(s) , (77) 



where £(s) is the Riemann zeta function. 

In order to compute the integrals in Eq. (|43p . it is convenient to perform an indefinite 
integration by parts and evaluate the boundary contributions only at the very end of the 
computation. 



/ 



t (2I) t (2I) r e t t (2I-l) t (2l) 

dt - t : = 7 + 21 dt --— = + t {2l) 



e* + e-* + 2 1 + e* J (1 + e*) e* 1 + e* 



dt * +et = -t (2l) Lio( - e*) + 2/ / dt t^ 2 '" 1 ) Li ( - e*) . (78) 

The procedure can be iterated thanks to the properties of the derivative of the polyloga- 
rithms. After 21 iterations we are left with 

21 fc-i 



/ dt - t ^ t 2 = -tW Li ( - e *) + E(-) fc+1 [ Ef (M " r)] Li fc ( - e *) . (79) 

fc=l r=0 



It is now easy to compute the boundary contributions. Since the integrand is an even 
function, it is enough to evaluate the integral in Eq. (|79p at t = and t — > —00 and then 
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doubling the result. By using the relations in Eq. (|77|) . one sees that at t = only the last 
term of the sum contributes, while all of the terms in Eq, (|79p vanish in the limit t — > — oo. 
Thus finally we find 

/ dt P t + P -t + o = — ( 2 ! Li 2 ;( - 1) = (2/)! (l - 2^ 2l A C(2i) ■ (80) 
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